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Abstract 

Purpose Multi-product processes are one source of 
multi-functionality causing widely discussed methodo¬ 
logical problems within life cycle assessment. A multi¬ 
functionality problem exists for comparative life cycle 
assessment (LCA) of multi-product processes with non¬ 
common products. This work develops a systematic 
workflow for fixing the multi-functionality problem 
caused by the non-common products. A novel technol¬ 
ogy for chlor-alkali electrolysis is analyzed and com¬ 
pared to the industrial standard technology to illustrate 
the approach and to benchmark the new technology's 
environmental impact. 

Methods A matrix-based workflow for comparative LCA of 
multi-product systems is presented. Products are distin¬ 
guished in main products and by-products based on the 
reason of process operation. We argue that only main prod¬ 
ucts form the reference flows of the compared multi-product 
systems. Fixing the multi-functionality problem follows di¬ 
rectly from the chosen reference flows. The framework 
suggests system expansion to fix the multi-functionality 
problem if non-common main products exist. Non¬ 
common by-products still cause a multi-functionality prob¬ 
lem. These by-products are systematically identified and the 
multi-functionality problem is fixed with avoided burden 
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and allocation. A case study applies the workflow for com¬ 
paring environmental impacts of the standard chlorine elec¬ 
trolysis to a novel process using oxygen-depolarized 
cathodes. Three scenarios are derived and evaluated. The 
assessed impact categories are cumulative energy demand, 
global warming potential, acidification potential, photo¬ 
chemical ozone creation potential, eutrophication potential, 
and human toxicity potential. 

Results and discussion The proposed workflow minimizes 
the methodological choices. The multi-functionality prob¬ 
lem is systematically fixed based on the distinction between 
the main products and by-products. Inconsistent solutions 
are prevented by rigorous identification of unequal by¬ 
products within the compared systems. Selecting avoided 
burden processes or allocation factors is the remaining am¬ 
biguous choice common to the standard methods. The case 
study demonstrates the applicability of the workflow to 
comparative LCA of multi-product systems. The case study 
results show lower environmental impacts for the novel 
electrolysis technology in all practically relevant scenarios 
and impact categories. 

Conclusions The framework for comparative LCA of multi¬ 
product systems with non-common products adds systemat¬ 
ic clarity to the general ISO standards. The approach 
reduces the subjective choices of LCA practitioners to the 
identification of reason of process operation. This reason is 
defined if the site-specific economic conditions are known. 
The matrix-based formulation allows identification of 
inconsistencies caused by multi-functionality. For the novel 
electrolysis technology, the results indicate significant po¬ 
tential for environmental impact reduction. 

Keywords Chlorine electrolysis • Comparative LCA • 
Industrial processes • Multi-functionality • Multi-product 
processes 
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1 Introduction 

The environmental impact of products, processes, and serv¬ 
ices can be analyzed using life cycle assessment (LCA). 
Environmental impacts are assessed using a functional unit 
as point of reference. The functional unit is fixed by LCA 
practitioners during the “goal and scope” phase of a LCA. 
For comparative LCA, the functional unit has to be carefully 
defined such that the compared alternatives fulfil an equal 
function. Severe methodological problems occur as soon as 
processes have more than one function (e.g., Azapagic and 
Clift 1999b, 2000). In many cases, the environmental 
impacts of these multi-functional processes have to be allo¬ 
cated to single functions. In LCA, this is referred to as the 
multi-functionality or allocation problem (e.g., Heijungs and 
Guinee 2007; Curran 2007; Finnveden et al. 2009). 

The ISO standard for LCA (ISO 2006a, b) provides a 
general code of practice for fixing the multi-functionality 
problem. The code of practice suggests a three-step proce¬ 
dure: The first step instructs to avoid allocation by two 
approaches: Dividing multi-functional processes into sub¬ 
processes or expanding the functional unit to include addi¬ 
tional functions. The latter procedure is called system ex¬ 
pansion. Among others, Ekvall (1999) and Weidema (2001) 
presented specific guidance on the application of system 
expansion. Still, allocation cannot always be avoided: 
Many multi-functional processes cannot be divided into 
independent sub-processes (e.g., chemical production pro¬ 
cesses). Expanding the functional unit may conflict with the 
goal of a study (e.g., environmental impacts of a single 
product). Therefore, the ISO standard suggests the applica¬ 
tion of suitable allocation methods in steps 2 and 3. 
According to step 2, such an allocation method should 
preferably represent an “underlying physical relationship” 
between the functions and the environmental impacts. 
Suitable relationships represent changes in inputs and out¬ 
puts caused by a change of a single function (Azapagic and 
Clift 1999a). As a last resort, the third step of the ISO 
standard suggests allocation based on other relationships 
between the functions, e.g., economic value. Application 
of economic allocation has been specified and discussed in 
several publications (e.g., Frischknecht 1998, 2000; Guinee 
et al. 2004). 

Still, the rather general formulations of the ISO standard 
(e.g., “wherever possible...”, “where allocation cannot be 
avoided...”) have been criticized (e.g., Heijungs and Guinee 
2007) because they leave a variety of methodological issues 
unsolved. In her literature review on the multi-functionality 
problem, Curran (2007) concluded that “the ISO standard 
should be expanded to provide more precise guidance in 
how to approach allocation.” In recent years, many authors 
give more specific guidance to fix the multi-functionality 
problem (e.g., Suh et al. 2010; Marvuglia et al. 2010). 


Moreover, several publications elaborate product-specific 
allocation procedures, e.g., for bio-fuels (Luo et al. 2009), 
agricultural products (Svanes et al. 2011), or combined heat 
and power production (Aldrich et al. 2011). 

This work aims at improving the applicability of the 
ISO standard for the case of comparative LCA of multi¬ 
product processes. In literature, multi-product processes 
are also referred to as multi-output processes (e.g., 
Curran 2007). Our motivation has its origin in a ques¬ 
tion raised in the chemical industry: Here, chlorine is 
mainly produced in the chlor-alkali electrolysis process, 
which jointly delivers gaseous chlorine, caustic soda in 
aqueous solution, and gaseous hydrogen as products 
(Schmittinger et al. 2002). Chlorine production is there¬ 
fore a popular example for a multi-product process in 
LCA literature (e.g., Boustead 1994; Frischknecht 1998; 
Weidema 2001; Schmidt 2009). A new electrolysis tech¬ 
nology has recently been developed to reduce 30 % of 
the electricity demand of the standard electrolysis tech¬ 
nology (Moussallem et al. 2008). This new technology 
uses oxygen-depolarized cathodes (ODC). The ODC 
process promises a lower electricity demand, but does 
no longer produce hydrogen. Comparing the ODC pro¬ 
cess to the existing chlor-alkali electrolysis process 
reveals a problem for industry and authorities: How 
are competing processes compared when they produce 
common (here chlorine and caustic soda) but also non¬ 
common products (here hydrogen)? Traditionally, LCA 
compares products. In LCA, the comparison of process¬ 
es is therefore reframed as comparison of multi-product 
systems. Comparing multi-product systems with non¬ 
common products raises the question whether the non¬ 
common products should be included in the reference 
flows of the functional unit. 

Even with a suitably defined functional unit, a multi¬ 
functionality problem remains for the case of non-common 
products. This can be rigorously shown by the analytical 
procedure developed by Heijungs and Frischknecht (1998), 
who studied multi-functionality from a mathematical point 
of view. They show that not all multi-product processes 
cause a multi-functionality problem depending on the cho¬ 
sen functional unit. These authors presented a rigorous 
analytical procedure to identify multi-functional systems 
that actually cause a multi-functionality problem. The pro¬ 
cedure shows that multi-functionality is an intrinsic problem 
for comparative LCA of multi-product processes with non¬ 
common products. 

To fix the multi-functionality problem in line with the 
ISO standard, we propose a systematic framework for com¬ 
parative LCA of multi-product processes with non-common 
products. The framework integrates both functional unit 
definition and fixing the multi-functionality problem. 
Herein, the procedure from Heijungs and Frischknecht 
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(1998) is further elaborated to not only identify a multi¬ 
functionality problem, but also guide the LCA practitioner 
to its solution in line with the ISO standard. 

The proposed framework is based on the matrix formu¬ 
lation of a LCA (Heijungs and Suh 2002). In Section 2, the 
fundamentals of this formulation are briefly summarized 
and the framework is developed. The framework is then 
applied in Section 3 to compare alternative chlorine elec¬ 
trolysis technologies that have non-common products. This 
case study is of major current interest to the chemical 
industry (Moussallem et al. 2009). The problem-specific 
aim of this work is to investigate the environmental impacts 
of the ODC process compared to the existing standard 
process. The results in Section 4 are therefore discussed 
from two perspectives: applicability of the proposed frame¬ 
work and interpretation of the case study results. 


2 Methods 

2.1 Problem formulation in matrix representation 

LCA of a linear process system can be represented by a set 
of matrices and vectors. This formulation is comprehensive¬ 
ly described by Heijungs and Suh (2002), whose notation is 
followed here. The next section briefly summarizes the 
basic fundamentals and highlights the scope of this work. 

For a given linear process system, the life cycle inventory 
(LCI) can be formulated as a set of linear equations defined 
by the technology matrix A, the scaling vector s, and the 
final demand vector f: 

A • s = f. (1) 

Environmental interventions of a system producing the final 
demand f are called LCI results. These results, summarized in a 
vector g, are computed from the intervention matrix B: 

g — B s. (2) 

The technology matrix A (economic flows) and the in¬ 
tervention matrix B (environmental flows) characterize a 
process system and are a result of data collection in LCA 
studies. The final demand vector f includes the reference 
flow of the functional unit. The functional unit is defined in 
the goal and scope phase of a LCA study. The reference 
flow may contain one or more economic flows. A scaling 
vector s is required to obtain LCI results from Eq. (2); this 
can be computed for an invertible and positive definite 
matrix A: 

s = A 1 f => g = B A 1 f. (3) 

In comparative LCA, a total number of r product system 
alternatives are analyzed. The matrices S, F, and G comprise 


sets of vectors s, f, and g, respectively, with columns repre¬ 
senting the £-th alternative: 

S = ( sj | ... | s k | ... | s r ); (4) 

F = ( f1 I ... | fk | ... | fr ); 

G= ( gi | ... | g k I ... | g r )• 

G = B S = B A 1 F. (5) 

Matrices A and B are assumed to contain all possibilities 
of implemented processes and are thus equal for all alterna¬ 
tive product systems. The environmental interventions in G 
are aggregated to potential environmental impacts by apply¬ 
ing common life cycle impact assessment (LC1A) methods. 
The LC1A methods are summarized in a characterization 
matrix Q. The LCI A results are finally collected in the 
impact matrix H: 

H = Q G = Q B A 1 F. (6) 

Multi-functionality is characterized by an over¬ 
determined system of linear equations. Such a multi¬ 
functional system has a non-square technology matrix A 
with more rows (i.e., economic flows) than columns (i.e., 
processes) and cannot be inverted. A multi-product process 
is consequently given by a process vector with at least two 
positive (i.e., output) entries that are seen as useful products 
(i.e., positive economic value). Computation of LCI results 
of multi-functional systems requires modifications of the 
technology matrix A and the final demand matrix F 
(Heijungs and Suh 2002). These modifications represent 
the procedures for fixing the multi-functionality problem 
and are commonly denoted system expansion and alloca¬ 
tion. However, system expansion can be applied in different 
variations. In the following, we refer to system expansion 
only if the functional unit is changed in the final demand 
matrix F. For clarity, we distinguish a modification that does 
not change the functional unit but accounts for multi¬ 
functionality with an avoided burden. We refer to this meth¬ 
od as avoided burden. 

Fixing the multi-functionality problem is not necessary 
for all non-square technology matrices A. Heijungs and 
Frischknecht (1998) presented a method that computes con¬ 
sistent LCI results from Eq. (2) with a non-square matrix A. 
The pseudo-inverse matrix A is used such that the absolute 
value of a discrepancy vector d is zero for some final 
demand vectors f: 

l|d|| = ||A • A + • f f|[ = 0. (7) 

Heijungs and Frischknecht (1998) concluded that a 
multi-functionality problem only exists if Eq. (7) is not true. 
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Recently, the use of the pseudo-inverse matrix A + was 
further elaborated to compute LCI results using different 
least-squares techniques (Marvuglia et al. 2010). These 
authors allowed the discrepancy vector d to be non-zero 
and computed LCI results without fixing the multi¬ 
functionality problem. 

In contrast, this work aims to fix the multi-functionality 
problem for comparative LCA of multi-product processes 
with non-common products. The discrepancy vector d is 
used to guide LCA practitioners through fixing the multi¬ 
functionality problem in line with the ISO standard. 
Following Fleijungs and Frischknecht (1998), the multi¬ 
functionality problem is considered to be fixed once the 
discrepancy vector d is zero. 

For a comparative LCA, the LCI results of multi¬ 
functional product systems can be computed from: 

G = B S = B A + F. (8) 

Here, the LCI results are only consistent if a discrepancy 
matrix D was a zero matrix: 

D mxr = (di | ... | d r ) = A • A + • F - F = 0" ,x '\ (9) 

In Eq. (9), the term A-A + -F can be considered as a final 
supply matrix F (cf. Heijungs and Suh 2002): 

F = A • A + • F. (10) 

If the final supply Fik is larger than the final demand Fik, 
there is an undesired surplus of the z-th product for the &-th 
system. These undesired surpluses are also identified by 
positive values for elements Dik = F,k — F^ in the discrep¬ 
ancy matrix D in Eq. (9). This has also been recognized by 
Marvuglia et al. (2010, p. 1028). 

2.2 Systematic framework for comparative LCA of multi¬ 
product processes 

A systematic framework for comparative LCA of multi¬ 
product systems is shown in Fig. 1. The procedure consists 
of two stages. In the first stage, the final demand vectors f k 
are defined. In the second stage, the multi-functionality 
problem is fixed. The discrepancy matrix D is used to 
identify the products causing a multi-functionality problem 
in both stages. 

In stage I, the first step of the suggested workflow 
ensures that the scope of a study is a comparative LCA of 
alternative multi-product systems. For the case of one multi¬ 
product LCA, the approach of Weidema (2001) can be used 
directly to fix the multi-functionality problem in line with 
the standard (ISO 2006a). The framework proposed in this 
work aims at studies where a total number of r multi-product 
systems are compared. Here, the reference flow typically 
consists of multiple product flows. A main focus is the case 


where multi-product systems have common and non¬ 
common products as it is often seen in the chemical indus¬ 
try. Here, the question arises whether the non-common 
products should be included in the reference flow of the 
final demand vectors. To answer this question, we suggest 
distinguishing main products and by-products for the com¬ 
pared multi-product processes. The main products are the 
site-specific reason of process operation. Remaining prod¬ 
ucts with a positive economic value are called by-products. 
The reason of process operation has been widely discussed 
(e.g., Oenning 1997). Similarly, Weidema (2001) uses the 
term “determining product” for the product that determines 
the production volume. The distinction between main prod¬ 
ucts and by-products varies depending on site-specific con¬ 
ditions and should always be defined in accordance with 
economic and technical boundary conditions. 

Based on the distinction of main and by-products, we 
suggest defining a system's function by the main product 
flows of the compared multi-product processes. 
Consequently, these main product flows should be used as 
reference flow in the final demand vectors f k . If the non¬ 
common products are main products, they should be includ¬ 
ed in the final demand vector. In this case, the number of 
main products in the vectors f k is unequal, i.e., the multi¬ 
product systems have unequal functions. We propose apply¬ 
ing system expansion for adjusting the final demand vectors 
f k to include equal numbers of main products. In this situ¬ 
ation, allocating environmental impacts only to the common 
main products is inappropriate from a practical perspective: 
If the non-common product is a main product, it needs to be 
produced by all compared multi-product systems. This prac¬ 
tical requirement is only captured by system expansion. 
However, system expansion requires additional processes 
for the production of the non-common main product. In 
practice, a total number of q competing processes may be 
suitable for system expansion. Accounting for all system 
expansion options increases the original number of r multi¬ 
product systems to (r+q~l) system alternatives, so that 

F=( fr I ••• I Wi). (11) 

If a process involved in system expansion is multi¬ 
functional, we suggest considering other than the main 
product(s) as by-products which do not change the final 
demand matrix F. 

It is possible that a process for system expansion does not 
exist. In that case, a comparative LCA of alternative multi¬ 
product systems is not reasonable because equal main prod¬ 
ucts cannot be delivered. Consequently, the scope of the 
study has to be reconsidered. 

While the compared systems should deliver the same 
main products after system expansion, comparability is not 
assured because the main products produced could be 
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Fig. 1 Systematic workflow 
for comparative LCA of multi¬ 
product processes 



supplied in different quantities. It is therefore useful to 
consider an adjusted technology matrix A MP . The matrix 


A mp contains only values for the main products (index MP) 
while all inputs and outputs other than the main products are 
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set to zero. A main product discrepancy matrix D MP can be 
calculated using the pseudo-inverse A + MP : 

Dmp = Amp • A^p • F — F. (12) 

The matrix D MP is used to verify that compared process 
systems deliver the main products in equal quantities. A 
non-zero discrepancy matrix D MP reveals unequal main 
product quantities for the alternative systems. Undesired 
amounts of a main product are identified by positive dis¬ 
crepancies (.Dmp);* in the discrepancy matrix D M p- In this 
case, the processes delivering the main products have to be 
analyzed for their flexibility in producing the main products. 
Flexible production of main products requires an adjustment 
of the process data in the technology matrix A and the 
intervention matrix B to reflect a change in process opera¬ 
tion (e.g., a combined heat and power plant with changing 
fuel demand and emissions for varying power-to-heat ra¬ 
tios). For the case of fixed production of main products 
(“joint production,” e.g., processes with fixed stoichiome¬ 
try), system expansion is again required to adjust the quan¬ 
tities of the main products. Stage I is concluded once a zero 
discrepancy matrix D MP is achieved and a set of comparable 
multi-product systems is established in the final demand 
matrix F. 

Stage II focuses on by-products of the multi-product 
systems. By-products may originate from the initial multi¬ 
product processes and from multi-functional processes 
added during system expansion. Non-common by-products 
can still result in an inconsistent solution of Eq. (8). To 
avoid an inconsistent solution, stage II continues with cal¬ 
culating the full discrepancy matrix D to check for a multi¬ 
functionality problem: 

D = A-A+-F-F=(d 1 |...| d r+? _i). (13) 

If the discrepancy matrix D is a zero matrix, there is no 
multi-functionality problem as concluded by Heijungs and 
Frischknecht (1998) and LCI results are computed consis¬ 
tently using Eq. (8). A non-zero discrepancy matrix D 
identifies now a multi-functionality problem caused by un¬ 
equal by-products. A positive discrepancy D ik marks a by¬ 
product (/-th row in D) and its corresponding multi-product 
system (Ar-th column in D) that needs to be modified to fix 
the multi-functionality problem. 

Before fixing the multi-functionality problem, we sug¬ 
gest identifying the by-product's site-specific use. Based on 
the site-specific use, the avoided burden method should be 
applied to account for the by-product surplus. The procedure 
represents technical reality most adequately (Weidema 
2001) and is mathematically described by Heijungs and 
Suh (2002, p. 41ff). 

But avoided burden processes can also be multi-functional. 
Therefore, recomputing the discrepancy matrix D is necessary 


using the adjusted matrices A burden and F burden . The recom¬ 
puted discrepancy matrix can be used to identify new 
unequal by-products introduced by the avoided burden 
process. Using multi-functional avoided burden processes 
may potentially cause an endless series of loops. An 
endless loop can be aborted by applying allocation. It is 
suggested to apply an allocation method once the eco¬ 
nomic or technical significance of newly introduced by¬ 
products is below a certain value. The significance can be 
checked by comparing LCI results using a suitable allo¬ 
cation factor with LCI results using an allocation factor 
that neglects the newly introduced by-product (i.e., the 
by-product has an allocation factor of zero). If the results 
differ only by a certain threshold value (e.g., max 1 % 
difference in relevant environmental interventions), then 
the loop can be aborted. 

In case an applicable avoided burden process does not 
exist, it is proposed to use an allocation method according to 
the ISO standard (ISO 2006b; Azapagic and Clift 1999a; 
Frischknecht 2000). The mathematical description of the 
allocation procedure is given by Heijungs and Suh (2002, 
p. 46ff). Allocation changes both technology matrix A and 
intervention matrix B by splitting columns according to an 
allocation factor. After allocation, a zero-discrepancy matrix 
D can be recalculated using the allocated technology matrix 
A a | lo and consistent LCI results can be computed from Eq. 
(8) using the allocated technology and intervention matrices 
A a iio an d Baiio- 

In the special case of all positive and equal discrepancies 
D ik in the /-th row, the /-th by-product is produced in equal 
amounts for all (r+q— 1) alternative systems. The surplus of 
this by-product can therefore be neglected in a comparative 
LCA. This is achieved by setting the respective process data 
Ay to zero for all processes. To summarize, the framework 
illustrated in Fig. 1 provides a systematic two-stage work¬ 
flow for the comparative LCA of multi-product systems 
with non-common products. 

3 Case study: technology change in chlorine industry 

3.1 Chlorine electrolysis technologies 

The so-called chlor-alkali electrolysis jointly produces gas¬ 
eous chlorine, gaseous hydrogen, and an aqueous solution 
of caustic soda from brine according to the following reac¬ 
tion: 

2 NaCl +2 H 2 0 -> Cl 2 +2 NaOH aq + H 2 . (14) 

The fundamentals of electrochemistry as well as the 
existing plant designs are not further elaborated in this work; 
the interested reader is referred to Schmittinger et al. (2002). 
The plant design considered as reference technology in this 
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study is the membrane process. The membrane process is 
labelled as the Best Available Technique for the chlor-alkali 
industry by the European Commission (2001) and common¬ 
ly used for new plants. 

A novel promising approach to reduce the electricity 
demand of the electrolysis process is using ODC (cf. 
Moussallem et al. 2008). The ODC process is based on a 
different overall chemical reaction: 

2 NaCl + H 2 0 +0.5 0 2 -> C l 2 + 2 NaOH aq . (15) 

This reaction requires 30 % less electricity than reaction 
(14) (cf. Moussallem et al. 2008). The difference is that 
hydrogen is not produced and pure oxygen is required as 
an input. Thus, membrane and ODC processes have only the 
products chlorine and caustic soda in common. Hydrogen is 
a non-common product only produced by the membrane 
process. The framework proposed in Section 2 is therefore 
applied to a comparative LCA of both technologies. 

3.2 System boundaries 

The system boundaries for the membrane process system are 
shown in Fig. 2a. The membrane process delivers gaseous 
chlorine and hydrogen as well as an aqueous solution of 
caustic soda of 32 mass percent. A concentration process is 
included to increase the caustic soda concentration up to the 
industrial standard concentration of 50 mass percent. 
Processes for electricity and heat generation are included 
to supply the energy demand of the electrolysis and the 
concentration process. Finally, the upstream impact of re¬ 
source supply, i.e., extraction of NaCl, is included. 
Figure 2b shows the ODC process system. The ODC pro¬ 
cess system additionally requires an air fractionation process 
for oxygen supply. The air fractionation process has nitro¬ 
gen as a by-product. A steam reforming process is included 
to account for hydrogen production of the membrane 
process. 

Auxiliary processes such as raw material preparation and 
product treatment are neglected in this work because data 
are not yet available for the new ODC process. For the same 
reason, manufacturing, installation, and deconstruction of 
the plants are also not considered here. Thus, the case study 
can be viewed as a simplified “cradle-to-gate” LCA for the 
products of both electrolysis technologies. The assumptions 
seem valid for a comparative LCA because the neglected 
processes and life stages are expected to have similar con¬ 
tributions for both membrane process and ODC process. 

3.3 Data sources, software, and impact assessment 

The LCA database ecoinvent (2011) contains LC1 data for 
the membrane process. No LCI data are available for the 


newly developed ODC process because it is not yet operated 
on an industrial scale. It is therefore necessary to determine 
the life cycle inventory data of the electrolysis with a theo¬ 
retical model. The electricity demand w e i of any electrolysis 
process depends on the voltage U and the specific Faraday 
constant/according to Faraday's law 

Wei = —,:■ (16) 

a 7 

Here, the electrolysis efficiency a accounts for line losses 
within the electrolysis process. Equation (16) is used to 
determine the electricity demands of both electrolysis tech¬ 
nologies depending on their voltages given in Table 1. The 
raw material demand is calculated using the stoichiometric 
reactions given in Eqs. (14) and (15). The heat demand of a 
concentration process can be calculated from Sattler (2001). 
The parameters used for modelling are summarized in 
Table 1. 

LCI data from the ecoinvent database for European aver¬ 
ages are used for modelling electricity and heat generation 
as well as NaCl extraction. Hydrogen production from 
steam reforming and oxygen supply from an air fraction¬ 
ation plant is modelled using LCI data from the GaBi 
database valid for Germany. The LCI data used for model¬ 
ling are summarized in Table 2. 

All scenarios are modelled using the GaBi software by 
PE International AG (2011). The impact categories cumula¬ 
tive energy demand (CED), global warming potential 
(GWPjoo), acidification potential (AP), photochemical oxi¬ 
dation potential (POCP), eutrophication potential (EP), and 
human toxicity potential (HTP) are evaluated using the 
CML characterization factors incorporated in the GaBi soft¬ 
ware tool. 

3.4 Application of the framework 

The problem-specific aim of this work is to investigate the 
environmental impacts the ODC process compared to the 
membrane process. This is done by applying the framework 
proposed in Section 2.2. First, the products of both technol¬ 
ogies have to be categorized as main and by-products. As 
indicated in Fig. 2, three possible main products (Cl 2 , 
NaOH, and H 2 ) are produced by the membrane process 
while the ODC process has only two possible main products 
(Cl 2 and NaOH). The actual main product definition 
depends on the actual site-specific conditions. In the follow¬ 
ing section, the framework is demonstrated for three practi¬ 
cally relevant scenarios. The scenarios differ in the chosen 
main products. The scenarios are distinguished by the su¬ 
perscript x of the corresponding final demand matrix FT 
The superscript x indicates the number of main products in 
the final demand vectors in Roman numbers. Membrane and 
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Fig. 2 Process flow sheets for 
a membrane process system 
and b ODC process system 


a) membrane process system 



b) ODC process system 
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ODC process systems are denoted by indices m and O for 
the final demand vectors f£, respectively. While the most 
important matrices are shown here, the Electronic 
Supplementary Material provides a detailed list of all matri¬ 
ces used in this work. 

The technology matrix A is shown in Table 3. It contains 
the processes shown in Fig. 2. Heat generation and concen¬ 
tration process for NaOH are equal for both electrolysis 
technologies and therefore not displayed to improve the 
readability of Table 3. 

The numerical values represent the stoichiometric mass 
relations according to Eqs. (14) and (15). The electricity 
demands are calculated from Eq. (16) using the parameters 
listed in Table 1. The matrix displays only the outputs of 
electricity generation, NaCI extraction, air fractionation, and 
steam refonner because aggregated LC1 data were used to 
model these processes (cf. Table 2). The environmental 
impact of the aggregated processes is included in the inter¬ 
vention matrix B (not displayed). 

Historically, the production of both chlorine and caustic 
soda has often been the reason of operation of chlor-alkali 
electrolysis plants. A scenario defining chlorine and caustic 
soda as main products for both processes appears natural. 
This leads to a final demand matrix F 11 shown in Table 3. 
Both alternative systems have an equal number of main 
products. The calculation of the main product discrepancy 
matrix D ^ p from the main product technology matrix 
Amp yields a zero matrix and thus concludes the main 


Table 1 Parameters used in modelling the electrolysis processes 


Parameter 

Name 

Unit 

Value 

U 

Membrane process 

V 

3.0 


voltage 



U 

ODC process voltage 

V 

2.0 

a 

Electrolysis efficiency 

- 

0.97 

f 

Faraday constant for CL 

(V-kg cl2 )/kWh 

1.3226 


product stage in the workflow (see Fig. 1), i.e., there are 
no non-common or unequal main products. 

Continuing with the by-product stage, the discrepancy 
matrix D 11 has two positive entries (Table 4). A surplus of 
the by-product nitrogen from the air fractionation process 
that supplies oxygen to the ODC process is indicated by 
Dq 4 > 0. The hydrogen surplus from the membrane pro¬ 
cess is identified by D I J l7 > 0. The surpluses are handled 
differently in this example. The steam reformer is used as an 
avoided burden process for hydrogen production. The nitro¬ 
gen surplus is treated by allocation because this example 
assumes that a nitrogen production process suitable for 
avoided burden does not exist. In fact, the LCI dataset for 
oxygen supply (cf. Table 2) is already allocated, so a manual 
allocation is not necessary. A recalculation of the discrep¬ 
ancy matrix D 11 using the adjusted technology and final 
demand matrices finally yields a zero matrix and allows a 
consistent computation of LCI results. 

A second scenario for defining the main products 
includes hydrogen as a third main product of the membrane 
process (cf. F IU in Table 3). This is especially relevant if 
hydrogen is required as a chemical and not available other¬ 
wise on the site of the chlor-alkali plant. The ODC process 
does not produce hydrogen. To establish comparability, the 
final demand vector of the ODC process system has to be 
expanded for hydrogen from the steam reformer (cf. f™ in 
Table 3). After system expansion, the main product discrep¬ 
ancy matrix D™ p is a zero matrix. The main products can 
be produced with equal quantities because the steam reformer 
process is scaled to the stoichiometrically fixed hydrogen pro¬ 
duction of the membrane process. The full discrepancy matrix 
D m in Table 4 yields only one positive value for Dq [ a . The 
nitrogen surplus is treated by allocation as described above in 
scenario F 11 before the LCI results can be computed. 

Electrolysis plants are traditionally built at the location of 
a chlorine demand. In contrast, caustic soda is often trans¬ 
ported to its place of use or stored if there is no demand. 
Therefore, the last scenario considers only chlorine as a 
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Table 2 Aggregated datasets 
used for upstream processes 
(cf. Fig. 2) 


RER: Europe; DE: Deutschland 
(Germany) 


Process 


Dataset 


Electricity generation 
NaCl extraction 
Eleat generation 
Air fractionation 
Steam refonner 


Ecoinvent: electricity, medium voltage, production mix, at grid (RER) 
Ecoinvent: sodium chloride, powder, at plant (RER) 

Ecoinvent: steam, for chemical processes, at plant (RER) 

GaBi-PE: oxygen, gaseous, from Linde process (DE) 

GaBi-PE: hydrogen, gaseous, from steam refonner (DE) 


main product. This example could also be assessed using the 
framework of Weidema (2001); however, it is examined 
here to demonstrate the flexibility of the suggested ap¬ 
proach. The final demand matrix F 1 of this scenario is 
shown in Table 3. A single main product for both alternative 
systems yields a zero main product discrepancy matrix 
D^jp. In the by-product stage, the full discrepancy matrix 
D 1 in Table 4 shows positive values for D I 0i and D' m 7 . 
These values indicate surpluses of nitrogen and hydrogen as 
already seen in scenario F 11 above. Additionally, D' m 6 > 0 
and Dq 9 > 0 identify the surplus of NaOH in both alterna¬ 
tive systems. The example in this work assumes an avoided 
burden process for NaOH does not exist. Therefore, alloca¬ 
tion is applied to both electrolysis processes and the air 
fractionation process before a zero discrepancy matrix D 1 
is found. In this example, we apply allocation based on the 
mass fraction of the products. However, the allocation 
method choice for chlor-alkali electrolysis is another 
discussion itself (cf. Boustead 1994) which is not further 
elaborated here. The Electronic Supplementary Material 
provides a brief discussion and results for alternative 
allocation methods. 


4 Results and discussion 

4.1 Results 

Table 5 summarizes the differences of the membrane 
(index m) and ODC (index O) process systems for six 
environmental impact categories h z . The results are pre¬ 
sented relative to the membrane process system illustrating the 
potential of the ODC process to reduce environmental 
impacts. 

For model validation, the results of the membrane 
process system were cross-checked with existing LCA 
data provided in the ecoinvent database. The results are 
in good agreement for the modelled life phases if the 
electricity demand used in the model of this work is 
adjusted to the ecoinvent dataset. The discrepancy of 
approx. 5 % between the data shown in this work and 
the ecoinvent dataset corresponds to the contribution 
from the following life cycle phases: raw material prep¬ 
aration, product treatment, transportation, and plant 
manufacturing processes. These life stages are neglected 
in this work (cf. Section 3.2) but included in the ecoin¬ 
vent dataset. 


Table 3 Technology matrix A and possible final demand matrices F 1 , F n , and F 111 for chlorine production scenarios 


Row Flow Unit A F 1 F 11 F 1 " 


Cl 2 Cl 2 , NaOH Cl 2 , NaOH, 

H 2 





Electricity 

generation 

NaCl 

extraction 

Air 

fractionation 

Membrane 

process 

ODC 

process 

Steam 

reformer 

fin 

f 1 

l o 

f n 

III 

f n 

l o 

.•III 

All 

r o 

1 

electricity 

kWh 

1 

0 

0 

-2.34 

-1.56 

0 

0 

0 

0 

0 

0 

0 

2 

NaCl 

kg 

0 

1 

0 

-1.65 

-1.65 

0 

0 

0 

0 

0 

0 

0 

3 

O2, air fractionation 

kg 

0 

0 

1 

0 

-0.23 

0 

0 

0 

0 

0 

0 

0 

4 

N2, air fractionation 

kg 

0 

0 

3.35 

0 

0 

0 

0 

0 

0 

0 

0 

0 

5 

CI2, membrane 

kg 

0 

0 

0 

1 

0 

0 

1 

0 

1 

0 

1 

0 

6 

NaOH mcmbrane 

kg 

0 

0 

0 

1.13 

0 

0 

0 

0 

1.13 

0 

1.13 

0 

7 

H2, membrane 

kg 

0 

0 

0 

0.03 

0 

0 

0 

0 

0 

0 

0.03 

0 

8 

CI2, ODC 

kg 

0 

0 

0 

0 

1 

0 

0 

1 

0 

1 

0 

1 

9 

NaOHoDC 

kg 

0 

0 

0 

0 

1.13 

0 

0 

0 

0 

1.13 

0 

1.13 

10 

f^ 2 , steam reformer 

kg 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0.03 


Heat generation and concentration process for NaOH are equal for both electrolysis technologies and therefore not displayed 
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Table 4 Full discrepancy ma¬ 
trices D 1 , D u , and D 111 for the 
different main product scenarios 


Row 

Flow 

Unit 

D 1 


D 11 


D 111 





<*L 

dg 

< 

dg 

d“ 

d m 

a o 

1 

Electricity 

MWh 

0 

0 

0 

0 

0 

0 

2 

NaCl 

kg 

0 

0 

0 

0 

0 

0 

3 

^2, air fractionation 

kg 

0 

-0.09 

0 

-0.21 

0 

-0.21 

4 

N2, air fractionation 

kg 

0 

0.03 

0 

0.06 

0 

0.06 

5 

CI2, membrane 

kg 

-0.56 

0 

0 

0 

0 

0 

6 

NaOH membrane 

kg 

0.50 

0 

0 

0 

0 

0 

7 

^2, membrane 

kg 

0.01 

0 

0.03 

0 

0 

0 

8 

CI2, ODC 

kg 

0 

-0.57 

0 

0 

0 

-0.02 

9 

NaOHoDC 

kg 

0 

0.49 

0 

-0.02 

0 

-0.02 

10 

^ 2 , steam reformer 

kg 

0 

0 

0 

-0.02 

0 

0 


4.2 Discussion of case study results 

The ODC process system has lower environmental impacts in 
all assessed impact categories for the scenarios F 1 , F 11 , and 
F m (cf. Table 5). The environmental impact reduction poten¬ 
tial of the ODC process system correlates strongly with the 
CED reduction potentials in all five assessed impact catego¬ 
ries: The coefficient of determination R 2 ranges from 0.90 
(HTP) to 0.99 (GWPioo) for the relative impact reduction of 
the ODC process system. The high level of correlation can be 
explained by the significant contribution of energy supply to 
the assessed impacts: The models in this work use the energy 
demand of the electrolysis processes and the concentration 
process exclusively to assess their environmental impact. The 
correlation with CED is not as strong for the HTP impact 
category. This impact category is more sensitive to the release 
of toxic substances which is minimized in large-scale energy 
supply systems. A comprehensive discussion of the generally 
observed correlation between CED and single-score LCIA 
indicators was recently presented by Huijbregts et al. (2010). 
The GWP reduction potential increases if an electricity mix 
with a higher GWP value is assumed. The steam reformer 
process does not depend on electricity generation; its relative 


contribution to the GWP thus decreases for electricity mixes 
with higher GWP values. 

The ODC process system leads to significantly larger 
reductions of relative impacts in the AP and POCP catego¬ 
ries. For these categories, the reduction potentials are larger 
than 15 % for both multiple main product scenarios F 11 and 
F m . In this case, the steam reformer process contributes 
only marginally to both impact categories. 

The presented results disregard direct emissions from the 
electrolysis plant that can have a significant influence par¬ 
ticularly in the HTP category. Neglecting this effect is 
justified in the comparison of membrane and ODC process¬ 
es at the present stage because no data are yet available for 
the ODC process. The development of the ODC process 
plant is aimed at reducing emissions compared to the exist¬ 
ing membrane process; thus, the result of the comparison 
should remain equal if these data are added. 

The contribution of auxiliary materials, NaCl preparation, 
transport and plant manufacturing, and recycling is also 
neglected in this work due to missing data for the ODC 
process. The environmental impacts of these items are 
expected to remain equal for both alternatives and shotdd 
therefore not strongly influence the result of the comparison. 


Table 5 Relative impact difference of membrane and ODC process system comparison using main product scenarios F 1 , F n , and F 111 



F 1 

Main product: CL; 
allocation: mass 

bio- 1 ^ 

bin, 

F 11 

Main products: CL and NaOH; 
burden: steam reformer 

h z n 0 -b” m 

h z,m 

pin 

Main products: CI 2 , NaOH, and Hn; 
expansion: steam reformer 

bfo-bii 

h“ 

CED 

-24 % 

-9 % 

-8 % 

GWP 100 

-22 % 

-3 % 

-3 % 

AP 

-24 % 

-20 % 

-19 % 

POCP 

-21 % 

-16 % 

-15 % 

EP 

-17 % 

-14 % 

-14 % 

HTP 

-10 % 

-11 % 

-11 % 
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The manufacturing and recycling of the plants might have an 
influence on the comparison because the ODC process plant is 
manufactured using different materials than a membrane 
plant. Manufacturing and recycling of the plant should be 
investigated in the future. However, it is expected that the 
manufacturing process of the plant does not significantly 
change the results: Existing data for the membrane process 
show contributions of only around 5 % in the assessed impact 
categories for this life stage (ecoinvent 2011). 

4.3 Discussion of systematic framework 

The case study demonstrates the advantage of the proposed 
framework: The methodological choices during a compara¬ 
tive LCA of multi-product systems are reduced to distin¬ 
guishing between the main products and by-products. This 
classification of main products and by-products can be seen 
as the question for the reason of process operation or the 
determining products (cf. Weidema 2001). Once such a 
classification is made, proper final demand vectors can be 
derived and existing methods for fixing the multi¬ 
functionality problem can be applied. These methods, name¬ 
ly system expansion, avoided burden, and allocation, natu¬ 
rally include further choices. Though the proposed 
framework cannot remove these choices themselves, it 
clearly indicates which method to use at a given point in 
the assessment. 

The benefit of the presented framework grows with in¬ 
creasing process integration (i.e., more by-product use) in 
future process systems. The integration of by-products often 
varies for different production sites. Therefore, the distinc¬ 
tion of main and by-products also depends on site-specific 
boundary conditions. The framework therefore represents a 
systematic step towards a site-specific application of LCA 
for comparative LCA. 

The application of the framework to the case study pro¬ 
vides a comprehensive view of the proposed methodology. 
Comparing electrolysis technologies may require choices 
for processes used for system expansion or avoided burden 
(e.g., steam reformer for H 2 production). Selecting from 
competing processes can be complicated and depends high¬ 
ly on the site-specific boundary conditions. The decision for 
a process used for both system expansion or avoided burden 
should thus be made in accordance with the scope of the 
study. For a site-specific assessment, the best available 
technology for the given site should be chosen. 

Multi-functional processes used as avoided burden may 
potentially cause an endless loop (cf. Fig. 1). As long as 
avoided burden processes of by-products are multi¬ 
functional themselves, choices for either allocation or new 
avoided burden processes are unavoidable. A rational ter¬ 
mination criterion has to be applied, e.g., as proposed at the 
end of Section 2.2. 


The proposed framework prefers system expansion and 
avoided burden over the application of allocation methods in 
line with the ISO standard. Even though all methods introduce 
uncertainty due to the number of possible choices (allocation 
factors, process selection), system expansion is closer to indus¬ 
trial reality. For comparative LCA, an alternative process will 
always be searched: an integrated production site would not 
replace a process as long as there are no alternative processes 
replacing possibly missing main products or by-products, even 
though an allocation-based study might advise them to do so. 

5 Conclusions and recommendations 

The presented work addresses the comparative LCA of 
multi-product systems, where alternative multi-product pro¬ 
cesses have common and non-common products. The pro¬ 
posed framework represents a systematic approach that links 
the functional unit definition to existing methods for fixing 
the multi-functionality problem. The final demand vector 
definition is reduced to the reason of process operation. This 
reason is usually given by economic boundary conditions 
and forms the basis for distinguishing main products and by¬ 
products. Hence, the final demand matrix is not an arbitrary 
choice, but rather a consequence of the scope of the study 
and site-specific economic conditions. 

Once the final demand matrix has been defined, the 
framework provides systematic guidance for fixing the 
multi-functionality problem. System expansion and avoided 
burden represent industrial reality better than applying allo¬ 
cation methods. The framework therefore proposes to either 
expand the functional unit (main product) or introduce an 
avoided burden (by-product). Allocation is applicable for 
by-products whenever reasonable avoided burdens are not 
available. The proposed framework does not allow an arbi¬ 
trary choice between avoided burden and allocation. The 
choices of the LCA practitioner are limited to choosing 
either suitable avoided burden processes or allocation fac¬ 
tors. The framework adds systematic clarity to the general 
formulation of the ISO standard focusing on comparative 
LCA of industrial multi-product systems. 

The presented work also points to future research in the 
field of multi-functionality in LCA. Possibly endless loops 
may occur whenever avoided burdens are multi-functional 
processes. The number of multi-functional processes is 
expected to increase in the future due to industrial develop¬ 
ment towards higher integration of material and energy flows. 
The presented methodology proposes a systematic approach 
to handle these loops. The proposed cut-off/allocation criteri¬ 
on should be further systematized in a mathematical approach 
to apply LCA to the design and optimization of future process 
systems. A promising approach in this direction was recently 
given by Marvuglia et al. (2010). 
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The results obtained in the case study show lower envi¬ 
ronmental impacts for the ODC process system in most 
scenarios and impact categories. To further strengthen the 
validity of these outcomes, the influence of direct plant 
emissions and plant manufacturing needs to be included in 
a next step as soon as data are available for the ODC 
process. The data in this work show already significant 
potential for the ODC technology to reduce the environmen¬ 
tal impacts of the chlor-alkali industry. 
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